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Abstract 

Rota's Umbral Calculus uses sequences of Sheffer polynomials to count certain 
combinatorial objects. We have developed a Maple package that implements Rota's 
Umbral Calculus and some of its generalizations. A Mathematica version of this 
package is being developed in parallel. 

Resume 

Le calcul ombral de Rota utilise des suites de polynomes de Sheffer pour enumerer 
certains objets combinatoires. Nous avons developpe une bibliotheque de fonctions 
Maple qui implementent cette theorie ainsi que ses generalisations. Une version 
Mathematica de ce package est developpee en parallele. 



AMS classification: 05A40. 

Hardware requirements: A computer equipped with Maple V.2 or V.3. 



1 Introduction 

Umbral calculus is the study of the analogies between various polynomial sequences and the 
powers sequence x n . For example, x n has many parallels with the lower factorial sequence 
(x) n = x(x — 1) • • • (x — n + 1): 

• The forward difference operator A : p(x) i— > p(x + 1) — p(x) plays a role with respect 
to (x) n analogous to that played by the derivative d with respect to x n . 
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Taylor's theorem is analogous to Newton's theorem. 



• The binomial theorem for (x + y) n is replaced by Vandermonde's identity for (x + y) n . 

Although Umbral Calculus dates back to the 18th century, it was only put on a rigorous 
foundation by Gian-Carlo Rota and his collaborators || [12[ in the 1970's. We now charac- 



terize each polynomial sequence under study by one or more polynomial operators (usually 
shift-invariantQ) associated with it. The duality between operators and polynomials is the 
key tool to deriving umbral calculus results. 

Umbral Calculus has many applications in enumerative combinatorics. The powers x n 
counts all functions from an n-element set to an x-element set, while the lower factorials 
(x) n count injections. Similarly, given any species of combinatorial structures (or quasi- 
species), let p n (x) be the number of functions from an n-element set to an x-element set 
enriched by this species. A function is enriched by associating a (weighted) structure with 
each of its fibers. The resulting sequence of polynomials {p n ) n& ]N * s sa id to be of binomial 
type since it obeys the "binomial" identity 

Pn(x + y) = J2 (^jPk(x)p n ~k(y)- 

For example, given the species of rooted forests, the enriched functions are called persistent 
functions and are enumerated by the Abel polynomials A n (x) = x(x + n) n_1 . Other 
applications include lattice path counting [[j], §, 

Our Maple package provides a number of different tools by which to enter operators. 
These operators can then be manipulated in many different ways. In particular, the poly- 
nomial sequences associated with them can be explicitly calculated. 

This package has already aided us in our research ||; we hope that it will help you too. 

We expect to release a Mathematica version of this package in the near future. 



2 Polynomial Operators 

Polynomial operators (shift-invariant or not) can be specified in several convenient man- 
ners: 

Explicitly by their action on polynomials. For example, the shift operator is defined 
< subs(x = x+a, p) | p > using the "angle-bracket" notation for functional oper- 
ators. (See ?operators [functional] and ?unapply for details.) Similarly, the 
Bernoulli operator p{x) \— > f£ +1 p(t)dt is defined <int(subs(x=t, p) , t=x. . (x+1)) 
I p I t> . 

As an analytic function of the derivative. By the expansion theorem (see |], Theo- 



rem 2, p. 185] or [^, Theorem 2, p. 691]), any shift- invariant operator can be ex- 
panded as a formal power series in d where d is a special reserved symbol representing 
the derivative. For example, the shift operator is defined exp(a*d). 



1 Shift-invariant operators commute with the shift operator E a :p(x) — > p(x ■ 



Abstractly as an unspecified function of d. For example, f (d) or f (d,x) in the case 
of a non-shift-invariant operator ||. 

Using the powseries package. If the coefficients of the formal power series given by 
the expansion theorem are all known, then use powcreate. For example, powseries 
[powcreate] (f (n) = a~n/n!);. 

As a series. If only finitely many terms are known, then use series. For example, 
series(l + a*d + a~2*d~2/2 + c*d~3, d, 3);. 

See ?linear for more information. 

Operators can be converted easily from one form to another with convert. A delta 
operator Q is a shift-invariant operator such that Q (x) is a non-zero constant. An ab- 
stract operator is assumed to be invertible unless indicated otherwise. (For example, Q : = 
d*f (d); or Q := f(d); f(0) := 0;.) 

Using the generalized expansion theorem, shift-invariant operators can be expanded 
with operatorExpansion as a formal power series in an arbitrary delta operator. 

> oe(d,delta,3) ; 

2 3 
exp(d) - 1 - 1/2 (exp(d) - 1) + 0((exp(d) - 1) ) 

Such operator expansions are practical for numerical calculations. For example, expanding 
the Bernoulli operator into powers of the forward difference operator delta yields the 
classical Newton-Cotes formula of numerical integration p. 186]. 

Our package also allows the expansion of linear operators which are not shift-invariant. 
Such expansions || express the operator as a formal power series in d whose coefficients are 
polynomials in x. For example, if Q is the operator Q:p(x) i— > f£ p(t)dt, then the expansion 
convert(Q, function, 5, x) or convert (Q , powseries ,x) of Q in terms of multiplication 
by x and the derivative d gives an elementary proof of Bourbaki's method of asymptotic 
integration |], Sections 3.5 and 3.6]. 

Operators can be applied to polynomials with dp by specifying the free variable of the 
polynomial: 

> dp(delta,x~3,x) ; 

2 

3 x + 3 x + 1 

In case the degree of the polynomial is not explicitly given, the program will calculate the 
Order most significant terms of the answer. (Order is a system variable whose default 
value is 6.) 



> dp(delta,x~n,x) ; 

(n-1) (n - 2) (n-3) 

n x + 1/2 n (n-1) x + 1/6 n (n-1) (n-2) x 



(n-4) 

+ 1/24 n (n-1) (n-2) (n-3) x 

(n-5) (n-6) 
+ 1/120 n (n-1) (n-2) (n-3) (n-4) x + 0(x ) 

3 Polynomial Sequences 

Given the necessary operators, the program can calculate polynomials of binomial type 
(bfo), Sheffer sequences (sfo), Steffensen sequences (steff), and cross-sequences (cseq) 
(see JL2|, Sections 5 an 8]). For example, the sequence of binomial type bf o(p(d), x, n) 
associated with a delta operator p(d) is defined by the conditions 

dp(p(d), bf o(p(d), x, n)) = n * bf o(p(d), x, n) for n>0 
bfo(p(d), 0,n) = for n>0 

bfo(p(d),x,0) = 1, 

or equivalently by its exponential generating function exp(q(x)t) where q is the composi- 
tional inverse of p. Note that q(t) is the generating function of the associated species. 

For example, the lower factorial (x) n is the basic sequence of binomial type for the 
forward difference operator A. 

> f actor(bf o(delta,x,4)) ; 

x (x - 1) (x - 2) (x - 3) 

If the degree is not explicitly given, then only the most significant terms will be computed. 

Several functions in the package do further operations on polynomial sequences. Ar- 
bitrary polynomials can be expressed in terms of such sequences (polynomialExpansion, 
shef f erExpansion, basicExpansion). For example, 

> p := randpoly(x); 

5 4 3 2 

p := 79 x + 56 x + 49 x + 63 x + 57 x - 59 

> be(delta,p,x) ; 

- 59 bfo(exp(d) - 1, x, 0) + 304 bfo(exp(d) - 1, x, 1) 

+ 1787 bfo(exp(d) - 1, x, 2) + 2360 bfo(exp(d) - 1, x, 3) 
+ 846 bfo(exp(d) - 1, x, 4) + 79 bfo(exp(d) - 1, x, 5) 

> "; 

5 4 3 2 

79 x + 56 x + 49 x + 63 x + 57 x - 59 

Connection constants can be determined between arbitrary polynomial sequences. For 
example, the Stirling numbers are given by cc(topseq(powerx, 5, x) , topseqQower , 
5, x) , x) where powerx(n,x) is x n and lower(n,x) is (x) n . Other features include 
umbral composition (uc), and umbral inversion (ui). 



4 Generalizations 



Several authors (eg. |T(J 14}| ) have generalized the umbral calculus by considering not 
only sequences of binomial type with generating function exp(g(x)t) but also those whose 
generating function is $(g(x)t) where = J2^=o^ n /N' an d W denotes the generalized 
factorial [n}\ = a(l)a(2) • • • a(n). Most of the functions in the umbral calculus package 
allow an optional argument a which is either left undefined, or defines the coefficients used 
by the "generalized derivative." Thus, 



> dp(d, x~3, x, proc(n) 1 



end) ; 
2 

x 



The following possible choices for a are predefined. 



Umbral Calculus 


dp(d,p,x,a) 


a 


Classical @, 12 
g-umbral calculus 10, 11 
Divided Difference 0, 13 

Hyperbolic || 


dp(x) 
dx 

p(qx) — p(x) 
(q — l)x 

p(x) — p(0) 

X 


classical (n) = n 
gaussian(n) = (q~n-l)/(q-l) 
divided(n) = 1 

hyperbolic (n) = 2*n*(2*n-l) 



15|| are supported. Most 



See ?genderiv for details. 

Generalizations of the umbral calculus to several variables 
functions included in the package have an alternate syntax for use in multivariate umbral 
calculi. In particular, d[i] represents the partial derivative with respect to the ith vari- 
able. Instead of a single delta operator, a collection of operators are required to define a 
sequence of binomial type. This generalization is completely compatible with the above 
generalization. See ?multi linear and ?moe for details. 

For further instructions consult the on-line help and examples provided in the package. 
For help, type ?key-word. An index of key-words is available via Tumoral. 

See for an extensive survey and bibliography of the umbral calculus. 
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